Terahertz inversionless lasing from a cavity-embedded two-dimensional electron gas 

in asymmetric quantum wells 
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Electric dipole transitions between different cavity polariton branches or between dressed atomic 
states with the same excitation number are strictly forbidden in centro-symmetric systems. For 
doped quantum wells in semiconductor microcavities, the strong coupling between an intersubband 
transition in the conduction band and a cavity mode produces two branches of intersubband cavity 
polaritons, whose normal-mode energy splitting is tunable and can be in the terahertz region. Here, 
we show that, by using asymmetric quantum wells, it is possible to have allowed dipolar transitions 
between different polaritonic branches, leading to the emission of terahertz photons. We present a 
quantum field theory for such a system and predict that high-efficiency, widely tunable terahertz 
lasing can be obtained. 



T3 

o 
o 



> 
m 
in 
o 

o 

(N 



X 



The resonant coupling between a photon mode and 
an electronic transition can produce hybrid light-matter 
eigenstates, such as the well-known dressed states in 
atomic cavity quantum electrodynamics ll[ or the po- 
lariton excitations in solid-state systems Q. In general, 
if one considers a doublet of dressed states, the electric 
dipole transition between these two states is strictly for- 
bidden in any centro-symmetric system. Analogously, 
in the case of cavity polaritons, a transition from the 
upper to the lower polariton branch cannot be accom- 
panied by any photon emission and can be only of non- 
radiative origin (e.g., via phonon emission Q). Recently 
though, it has been shown that, in exciton-polariton sys- 
tems, the application of an electric field produces hy- 
bridisation of internal motion exciton states with differ- 
ent parity, thus breaking such a selection rule This 
enables transitions between different exciton-polariton 
states and in principle paves the way to integrated tera- 
hertz (THz) sources where the excitation is provided by a 
near-infrared pump with potentially significant quantum 
efficiencies (up to 1.5%). 

Indeed, sources emitting in the THz region of the elec- 
tromagnetic spectrum are a subject of intense investiga- 
tions both for fundamental physics and applications. In 
fact the THz region represents a technological gap that 
can be covered only with difficulty, either with electronic 
devices (at lower frequencies), or with lasers (at higher 
frequencies). A great research effort has been accom- 
lished to extend semiconductor quantum cascade lasers 
to the THz region in such devices transitions oc- 
curs from quantized subbands in the conduction band of 
suitably designed multiple quantum well structures. Al- 
though THz quantum cascade lasers are developing 
fundamental limitations occur because the radiative life- 
times of the excited electronic states are long compared 
to the intrinsic fast non-radiative recombination chan- 
nels that arise from efficient electron-phonon scattering 
(at least in the case of incoherent electronic excitations 



In the case of intersubband transitions, it is possi- 
ble to reach the strong light-matter coupling regime by 
embedding a structure containing mu ltip le doped quan- 
tum wells in microcavity resonators Ill4l4|. leading to 
the creation of so-called intersubband cavity polariton 
modes. In such a system, it is possible to control in-situ 
the value of the polariton splitting and it is pos- 

sible to engineer the intersubband transitions in a very 
flexible way [ljj ]. A fundamental unexplored question is 
whether in such a class of intersubband systems it would 
be possible to have efficient radiative transitions between 
polariton branches. 

In this letter, we show that breaking the wavefunc- 
tion symmetry using asymmetric quantum wells enables 
radiative transitions between different intersubband po- 
lariton branches, leading to intrinsically efficient, widely 
tunable THz lasing without electronic population inver- 
sion. We develop a quantum theory describing such pro- 
cess and determine the overall efficiency of the THz emis- 
sion, in a way that takes into account the nonbosonicity 
of intersubband excitations 0, El • We have calculated, 
both analytically and numerically, the intrinsic quantum 
efficiency for the considered process showing that, for 
realistic parameters, unprecedented quantum efficiencies 
are achievable. 

Let us start by considering a quantum Hamiltonian de- 
scribing two parallel semiconductor conduction subbands 
coupled to a quasi-resonant microcavity photon mode, in 
the rotating wave approximation (RWA), namely 



H = 
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where k and a q are, respectively, the creation oper- 
ators for an electron in the quantum well conduction 
subband [i with in-plane wave vector k and a micro- 
cavity photon with in-plane wave vector q. Their re- 
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FIG. 1: Typical energy dispersion of the lower and upper in- 
tersubband cavity polariton branches (black solid lines), and 
of the THz mode of the metallic cavity (red dashed lines). 
All the allowed final states for the radiative process under ex- 
amination, in which an upper polariton scatters into a lower 
polariton plus a THz photon, are marked in red on the 3D 
image. For clarity the lightcone of the THz mode is also 
shown originating from the upper polariton mode. The pro- 
cess highlighted by red dots corresponds to the forward scat- 
tering process considered in the text, such that q and q' are 
parallel. 



spective energies are huj^(k) and hw-p^q), with hi02(k) = 
fvjj\(k) + fvjj%2. The coefficient hx(q) quantifies the light- 
matter coupling and it is proportional to the intersub- 
band dipole z\i — f dz zijji(z)-ijj2(z), where ip^ is the 
wave function along z, normal to the quantum well plane, 
of an electron in subband fj,. If we consider an inter- 
subband energy Swi2 in the mid-infrared region, the 
Coulomb interaction 



13l . l2fj, |2l| can be neglected to a 



first approximation. All the interactions are spin con- 
serving and are between states in the same quantum well, 
so for simplicity we omit the spin and quantum well in- 
dices. Due to the intersubband transition selection rules, 
only the Transverse Magnetic (TM) cavity mode couples 
to electrons. Introducing the (approximately bosonic) 
bright intersubband operator 6 q [12j, 



K = 



i 



k 



2,k+q c l,kj 



(2) 



the intersubband dynamics from Eq. ([T]) can be modeled 
by the effective Hamiltonian 

H' = h ^2 w ph(g) a q a q + wi 2 6 q 6 q + £l(g)(a q 6 q + & q a q ), 
q 

( 3 ) 

je{LP,UP},q 

where N = riQw-^2DEG<S' is the total number of electrons 



in the structure (uqw is the number of quantum wells, 
iV2DEG is the density of the two-dimensional electron gas 
and S is the transverse surface) and fl(q) = VNx(q) 
is the vacuum Rabi frequency. The Hamiltonian in Eq. 
([3]) can be diagonalized in terms of the normal-mode po- 
lariton operators pt = Xj^a^ + j/j,g6 q , where Xj >q and 
Dj. q are the (real) Hopfield coefficients that denote the 
light and matter fractions of the polaritonic excitations 
respectively. In Fig. [TJ a sketch of the typical energy- 
momentum dispersion for the lower and upper polariton 
branches is shown. 

For intersubband cavity polaritons with energies in the 
mid-infrared, transitions between the upper and lower 
polariton branch are typically in the THz region of the 
spectrum. The cavity also contains photon modes at 
these THz frequencies, and here we consider the process, 
depicted in Fig. [TJ where a polariton is first excited into 
the upper branch with a pump beam, and then scatters 
to the lower branch by emitting a THz photon into the 
cavity. Such a coupling to the the low energy THz modes 
is described by the interaction term: 



-ffint = ef 




(a q e ifq 



*), (4) 



where L cav is the microcavity thickness, f the electron 
position operator and cv q is the creation operator of the 
low energy THz cavity photon with energy Swthz (<?) ■ For 
definiteness, in this letter we will consider the case of a 
metallic cavity, in which cv q is the TEM mode with lin- 
ear dispersion [14j . while a q belongs to an excited cavity 
branch, but other configurations are possible. 

The matrix element describing the process illustrated 
in Fig. [1] in which a polariton jumps from the upper to 
the lower branch by emitting a THz photon, reads 



(G|a p p L p,q' .ffint Pup JG) = y VP , q y hP , q ,eAzJ ^™ z ^ 

x%-q'-p), (5) 

where |G) = rifc</c F c i kl^) is the Fermi ground state 
of the electronic system, |0) is the vacuum state and 
Az = J dzz [\ip-2.{z)\ 2 —\^pi{z)\ 2 ] is the interbranch dipole, 
representing the distance between the mean electron po- 
sitions in the two electronic subbands. As we have antic- 
ipated, this scattering process occurs only in asymmetric 
systems, where Az ^ 0. 

It is apparent from Fig. [T]that, for a given initial up- 
per polariton mode, there is a continuum of final modes 
conserving both energy and momentum, given by the in- 
tersection between the lightcone of the THz mode irra- 
diating from the initial upper polariton mode and the 
lower polariton branch. From Eq. ([5]) we see that the 
final mode with the larger value of q' , that is the one 
with the highest matter fraction yhp, q ', is generally the 
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one most likely to win in a mode selection competition in 
a THz cavity. For definiteness then, in the following, we 
will confine our attention to just the forward scattering 
process shown in Fig. [TJ in which q and q' are parallel. 
Now we are concerned only with three modes; the up- 
per polariton, the lower polariton and the THz photon. 
Their in-plane wave vectors are fully determined by the 
argument above, so we can simplify the notation omit- 
ting the wave vectors for the quantities relative to these 
modes. 

The spontaneous emission matrix element in Eq. ([5]) 
is not enough on its own to let us calculate the THz 
emission in the stimulated regime; instead we need the 
full many body matrix element describing a transition 
from a state with ]Vtjp upper polaritons, ATlp lower po- 
laritons and -/Vthz THz photons to a state in which 
an upper polariton has become a lower one, emitting 
a THz photon. Were intersubband excitations elemen- 
tary bosons, the transition rate would be exactly equal to 
the spontaneous rate times the bosonic enhancement rate 
N\jp(Nlp + 1)(JVthz + !)■ However, intersubband exci- 
tations behave as composite bosons (l9j . implying that 
the bosonic enhancement factor is renormalized by the 
bosonicity factor -B^p < 1 which depends on the popu- 
lations of the initial state Q . For normalized excitation 
densities jVup j ^ jVLF -C 1, this coefficient can be well ap- 



proximated by B^t hF — 1 — C: 
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where Clp 



and Cup , each of the order of one, depend on the Hopfield 
coefficients of the polariton modes. Notice that a state 
with polaritonic populations N\jp and TVlp corresponds 
to a linear superposition of states which has, at most, 
iVup + A^lp electrons in the excited subband. 

Denoting the linewidths of the different modes as Typ, 
Tlp and Tthz, and assuming, for simplicity, that the 
lower polariton has a Lorentzian lineshape, the sponta- 
neous scattering rate 5 S for the considered interbranch 
dipole transition can be approximated via the Fermi 
golden rule as 



^THz(Aze^upy L p) s 
2her > e r SL c&v 



r LP 5(w UP - W T Hz ^ w) 



In the limit -B^ LP ~ 1, we can solve Eq. ([8j for its 
steady state and compute the quantum efficiency as 



( W - WL p) 2 + (^) 2 

where we have assumed, as is usually the case experimen- 
tally, that the width of the lower polariton mode is much 
larger than that of the THz photon mode (r L p 3> r THz )- 
For the resonant, stimulated process, this becomes 



(6) 



2B™u Tm (Azey UP y LP y 

heoe r SL caiV rLp 
1) - (N VP + l)N LP N TRz ] 



(7) 



E(N V p,N L p,N THz ) 
x[iV UP (7V LP + l)(iVTHz 

If we denote the pumping rate into the upper polariton 
mode as P, we can thus write the following rate equation 
for the populations in the three modes 

a\jp = -r U p^ UP -s(7v UP ,iv L p,^ T Hz) + P 

NhP - -T L pN L p + E(N V p,N L p,N Tliz ) (8) 
A^THz = -rTHzA r THz + S(Ar U p,A^ L p,7V TH z). 



c UP^ THz ^ THz 



max[l 



,0] 



1 



r L p 



(9) 



where the Hopfield coefficient x\jp accounts for the fact 
that the pump beam couples only to the photon part of 
the upper polariton. From Eq. §§§ it is possible to verify 
that, as we anticipated, extremely high quantum efficien- 
cies can be achieved going toward higher pump powers. 
We can also write the condition for stimulated THz emis- 
sion to be possible as S(A r yp, 0, 0) > Tthz- From Eq. ([7]) 
the threshold density for stimulated emission, neglecting 

b °n up - is thcn 



N\jp _ Tthz 7leoe r £ cav rLp 
S Wth z 2(Azey V py- L p) 2 ' 



(10) 
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FIG. 2: Black solid line: quantum efficiency r\ as a function 
of the applied pump power density calculated solving numer- 
ically Eq. (J3J, including thus nonbosonicity effects. Black 
dashed line: same quantity obtained using Eq. ((9|. Red 
dash-dotted line: total excitation fraction UP jv LF • Inset: 
wavefunctions for the first two conduction subbands for an 
asymmetric GaAs quantum well structure. Due to the quan- 
tum well asymmetry, Az — (ip2\z\ip2) — {"^\\z\ij)i) 7^ 0. All 
the parameters are given in the text. 



Using this formalism to predict quantitative conver- 
sion efficiencies requires us to consider a specific exam- 
ple structure. In general, a compromise is needed to 
get a large enough Z\% to achieve the strong-coupling 
required to create a polariton frequency splitting in the 
THz range and a large enough Az to have an interbranch 
transition dipole. As an example we consider the simple 
asymmetric stepped GaAs quantum well of Fig. [2] with 
infinite barriers. It is characterized by only three pa- 
rameters, the overall width, Lqw, the height, V, and 
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the position, D, of the potential step. The wavefunc- 
tions shown in Fig. [2] correspond to £qw = 25 nm, 
D = 0.4L QW and V = 138 meV; this gives an ISBT en- 
ergy huji2 ~ 100 meV and Az ~ z 12 ~ 0.1Lq W . These 
values are not optimal, but similar ones are obtainable 
in a fairly large sector of the parameter space and we 
thus expect them to be easily achievable also in more 
realistic geometries. We model a structure composed of 
riQw — 40 such GaAs QWs each doped at an electron 
density N 2 deg — 10 12 cm~ 2 , that gives a resonant vac- 
uum Rabi frequency Q(q rcs ) ~ O.IW12 [where g lcs is the 
resonant wave vector such that 0J p h(<Zres) = ^12]- As- 
suming j/up — 0.6 and j/lp — 0.9 (the values of the 
process represented in Fig. Q]), a quality factor for the 
THz mode = 100 [22| and polaritonic linewidths 

Tup = T LP = 4 meV 0, Eq. (JTUJ) gives ^ ~ 0.05, 
that is to say that the onset of stimulated THz emis- 
sion occurs when only ~ 5% of the electrons are in the 
excited subband. Such low excitation threshold density 
justify a posteriori the approximation of neglecting the 
nonbosonicity factor in Eq. (fT0|) . 



In Fig. [3] we explore the versatility of this system 
by studying the way the forward-scattered emitted THz 
photon frequency wthz depends on the pump frequency 
wup, for different values of the light-matter coupling 
^(<?res)- At lower values of Q(q TCS ) the mechanisn gives a 
frequency down-conversion up to a factor 10, while stay- 
ing in the strong coupling regime and maintaining large 
conversion efficiencies. Since Q,(q ICS ) can be modulate d by 
electronically changing the electron gas density HI 17j . 
this offers a way of realizing widely tunable THz emitters. 

In conclusion we have shown how, through the use of 
asymmetric QWs in intersubband polariton systems, it 
is possible to obtain stimulated THz emission character- 
ized both by an extremely large quantum efficiency and 
a remarkable frequency tunability. 
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FIG. 3: Frequency of the emitted THz radiation wthz, consid- 
ering forward scattering, as a function of the pump frequency 
cjup- The different lines correspond to = 0.05 (solid line), 
0.1 (dashed line) and 0.15 (dash-dotted line). The red circle 
marks the parameters considered for the numerical applica- 
tion in the text. 

This offers the possibility of demonstrating a THz in- 
versionless laser [2^, [24[. In Fig. [5] we see the quantum 
efficiency as a function of the pump power obtained by 
numerically solving Eq. © (solid black line) and using 
Eq. ([5]), that is, neglecting nonbosonicity effects (black 
dashed line). In the same Figure we also plot the total 
excitation density Nvp + Nl - P as a function of the pump 
power (red dash-dotted line). We see that, with exper- 
imentally achievable pump powers [2^, a quantum con- 
version efficiency in excess of 20% is possible, still with 
the system being a long way from population inversion. 
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